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Abstract 

The energy production through thermo-nuclear fusion requires the confine- 
I ment of the plasma into a bounded domain. In most of the cases, such config- 



urations are obtained by using strong magnetic fields. Several models exist for 
describing the evolution of a strongly magnetized plasma, i.e., guiding-center ap- 



^ ■ proximation, finite Larmor radius regime, etc. The topic of this paper concerns 

!>■: 

QQ . a different approach leading to plasma confinement. More exactly we are inter- 

■ ested in mathematical models with fast oscillating magnetic fields. We provide 

, rigorous derivations for this kind of models and analyze their properties. 

(N 
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1 Introduction 

Motivated by the energy production through thermo-nuclear fusion, many research 

programs concern plasma confinement models. It is well known that good confinement 

properties are obtained under strong magnetic fields 5^ = 0{l/e) with e > a small 
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parameter. Using the kinetic description and neglecting the particle collisions lead to 
the Vlasov equation 

dtf + — ■ Vxf + e (E'{t, x) + — A B'{t, x)) ■ Vpf = 0, {t, x, p) G M+ x x (1) 
with the initial condition 

r(0,x,p) = r(x,p), (x,p)gM3xR3. (2) 

Here = f'^{t^x^p) > is the distribution function of the particles in the position- 
momentum phase space G R'^ x M^, m is the particle mass, e is the particle charge 
and (i?^, B'^) stands for the electro-magnetic field. 

Standard configurations ensuring confinement are those obtained by applying strong 
magnetic fields. For example, assuming that the electric field derives from a given 
potential E = —'Vx4> and the magnetic field is stationary, divergence free 

B(x) 

B%x) = -^b{x), div^iBb) =0, < £ « 1 

for some scalar positive function B{x) and some field of unitary vectors b{x), lead to 
the Vlasov equation 

dtf + - ■ V.r + (eE{t, x) + A b{x)] ■ Vpf = 0, Uc{x) = (3) 

m \ e J m 

whose limit as e \ is known as the guiding-center approximation. The particles 

rotate around the magnetic lines and the radius of this circular motion, which is called 

the Larmor radius pL, is proportional to the inverse of the magnetic field. Therefore 

when the magnetic field is strong, the typical Larmor radius vanishes and the particles 

remain confined along the magnetic lines. But the frequency of these rotations, which is 

called the cyclotronic frequency, is proportional with the magnetic field. Consequently, 

high magnetic fields introduce small time scales, since the cyclotronic period is much 

smaller than the observation time. Clearly, the transport equation ([3D involves multiple 

scales: fast motion around the magnetic lines driven by the Laplace force in ^^^^^(p A 

b) ■ Vp and slow motion corresponding to the advection ^ ■ V^,. + eE -Vp. 

From the numerical point of view, the efficient resolution of requires multiple 

scale analysis [3] or homogenization techniques. It is also possible to appeal to La- 

grangian and Hamiltonian methods For a unified treatment of the main physical 



ideas and theoretical methods that have emerged on magnetic plasma confinement we 
refer to | T3] . 

The guiding-center approximation for the Vlasov-Maxwell system was studied in 
[5] by the modulated energy method. The case of three dimensional general magnetic 
shapes ([3D has been studied in [7], using a general method, based on ergodicity, in- 
troduced in [6J. It was proved in [7] that the limit density / = lime\^o satisfies the 
Vlasov equation 

dtf + b(x) (g) b(x)— ■ Vxf + (eb(x) (g) b(x)E + uj(x,p) p) ■ Vpf = 
m 

where for any {x, p) with pAb{x) ^ the symbol p stands for the orthogonal momentum 
to p, contained in the plane determined by b{x) and p, and such that its coordinate 
along b{x) is positive, that means 

p=\pA b{x) I b{x) - {p ■ b{x)) 

and the frequency uj{x,p) is given by 

, , \pAb(x)\ , (p ■ b(x)) f db , , . p \ ,/ x 

^ix,p)= ^ ^' div,.6- ^ I b{x)-—^), pAb{x)j^O. 
2m m \ox \pAb[x)\/ 

The analysis of the Vlasov or Vlasov-Poisson equations with large external magnetic 

field have been carried out in [10], [12], [8], [TT], [13]. The numerical approximation 

of the gyrokinetic models has been performed in [H] using semi-Lagrangian schemes. 

Other methods are based on the water bag representation of the distribution function 



Notice that configurations with large magnetic field amplitude require huge energy 
since the magnetic energy is quadratic with respect to = B/e. 
We investigate here models with fast oscillating magnetic fields 

B'{t, x) = 9{t/£)B{x)b{x), < e « 1 (4) 

where 6 = 6{s) is a given T periodic profile of class C^. The magnetic energy dissipated 
in this case is much lower than for the guiding-center approximation and remains of 
order of \B\'^. Therefore such models will be more interesting for real life applications, 
provided they still have good confinement properties. At this stage we neglect the 
gradient and curvature effects of the magnetic field, assuming that 

B' = B'it) = e{t/e){0,0,B) 



for some constant B > 0. The general model including gradient and curvature effects 
will be discussed in Section [71 The vector potential corresponding to this magnetic 
field, i.e., satisfying = cur^A^, div^.^!^ = is given by 

A%t, x) = -— e{t/e)^x, ^x = (x2, -xi, 0). 

Decomposing the electric field into gradient and rotational parts 

= -V^(j) + cml^^' 

we deduce, by Faraday's law dtB"^ + cur^i?^ = that 

cuil^idtA" + cmlr^ip^) = 0, div^((9iv4^ + curU?/;^) = 

and therefore the electric field induced by the time fluctuations of the magnetic field is 

cmhr = -dtA' = ^e'{t/e)^x. 

The Vlasov equation becomes, with the notations -^p = {p2, —pi, 0) and E = —Vx(p 

dtf + £ ■ + {eE{t, x) + '^e\t/e)^x + e{t/e)^p) ■ VJ' = 0. (5) 

Here E = — Va;0 is a given irrotational electric field or can be determined eventually 
by solving the Poisson equation 

div,E(t) = -A,0(t) = - ( / f{t,x,p)dp-no{x)\, tGM+, xGMl (6) 

The concentration no(x) corresponds to a neutralizing background of charged particles 
of opposite sign and Eq is the electric permittivity of the vacuum. 

Our paper is organized as follows. The main results are presented in Section |2J 
In Section [3] we introduce the mathematical tools that we need for our analysis. It 
mainly concerns the average operator with respect to characteristic flows. We discuss 
its main properties as range characterization, Poincare and Sobolev inequalities. Some 
commutation results are established in Section HJ Section Ois devoted to the derivation 
of the limit model, which follows in a natural way by appealing to the average operator 
introduced before. We establish the conservation of the total energy and justify the 
confinement properties for such a model. The asymptotic behaviour towards this limit 
model is analyzed in Section |6l The general three dimensional setting is investigated 
in the last section. 
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2 Presentation of the model and main results 



The Vlasov equation ([5]) reduces to the characteristic system 

^ = ^ = eE{t, X^{t)) + ^9 \t/e) ^X%t) + coAt/e) ^P%t). (7) 

It is convenient to introduce the fast variable s = t/e and the standard ansatz 

X'{t) = X\t,t/e) +eX\t,t/e) + P'{t) = P\t,t/e) + eP\t,t/e) + ... 

(here all dependences with respect to the fast variable s are supposed T periodic, as 
the profile 6 = 6{s)) leading to 

dtX'' + -dsX' + sidtX' + -dsX') + ... = — + e— + ... 

e e mm 



and 



dtP' + -dsP'' + e{dtP' + -d,P') + ... = eE{t,X'' + eX' + .. 
e e 



+ '^e'{t/e)^{X' + eX' + 
uj,e{t/e)^{P^ + eP^ + ...). 



At least formally one gets the equations 



9,X° = 0, 9,P° = ^^'(s)^X0 (8) 

at the lowest order and 

5^X0 + dsX^ = — , dtP"" + dsP^ = eE{t, X°) + "^e'is) ^X' + io^s) ^P' (9) 
m 2 

at the next order e^. It follows that the quantities 



depend only on t. Therefore, in order to obtain the characteristic equations satisfied 
by the leading order terms (X°,P°) we write the equations in terms of {X^,Q^) 
and eliminate {X^, P^) by averaging with respect to the fast variable s over one period. 
The first equation in ([9]) becomes 

m 2 
5 



and therefore averaging with respect to s yields 



dt m 

Similarly, the second equation in implies 



+ ^(0) ^X°(t), {0) = ^l 0{s)ds. (10) 



+ u;As)HQ' + ^e{s)^X'} 
= eE{t,X'^) + ^Ue^X'} 



2 ' ' [m 2 

+ uj,e{s)^Q'+'^e\s)^^x\ 



Finally one gets 



and therefore, averaging with respect to s yields 

^ = eEit, X\t)) + I {9) ^Q^t) + ^ {9') ^ ^X\t), (9^) = 9\s) ds. 

(11) 

We associate to the characteristic equations ( ITU]) . (ITT]) the transport equation 

dtg' + + Y {9) ^x) ■ V../ + (eE(t, x) + | {9) + ^ {9') ^^x) ■ W,g' = 0. 

(12) 

Since we have 

{X%t),P%t)) ^ {X\t,t/e),P%t,t/e)) = [x%t),Q\t) + ^9{t/e) ^X%t)) 

and assuming that X'^{0) = x, P'^(O) = p we can write 

f{t,X\t,t/e),P\t,t/e)) ^ r{t,X%t),P%t)) = r{x,p) 

g\t,X\t,t/e),P\t,t/e) - -^9{t/e) ^X\t,t/e)) = g\t, X\t),Q\t)) = g'^ix,q) 
and therefore we can expect that 

r{t,x,p) ^ g^{t,x,q = p ^9{t/e) ^x), as e \ 

provided that the initial conditions are well prepared. Introducing the density /°(t, s, x, p) 
g^{t,x,p — mijJc9{s)/2 -^x) in the phase space {s,x,p) we deduce that f^(t,x,p) ~ 
/"(t, t/e, x,p) as e \ 0. 
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Theorem 2.1 Assume that E e Aoc(^+; L°°(^^)), f e L'^{R^ x R^). For any e > 

let G L°°(]R_|_; x M^)) &e a weaA; solution of Then there is a sequence {6n)n 

converging to zero such that {f'^")n two-scale converges towards a weak solution of 

dtf+{^-'^{e{s)-{e))^x).vj' (13) 

+ {eE{t, X) + ^(^(.) + {9)) + ^{{e^) - e\s)) ^^x) ■ V,f = 

f\0,s,x,p) = f'^ (x,p-'^{e{s) - em ^x) ekeiT. (14) 

Consequently we have to study the confinement properties of the limit model (fT2i) 
(or (|T3l) ). Indeed, such models lead to confinement. For convincing ourselves let us 
consider a particular case, that of vanishing electric potential = 0. The characteristic 
system for (X°, Q^) becomes 

implying that 

- - W') ^^X%t) = {6) (15) 

dX^ 

Multiplying by we obtain the conservation 

dt 

\^m-{9f)l\^^X^^' 
If 6 is not a constant profile {i.e., the magnetic field oscillates in time), then 

%m-{of) = ^{{e-{e)r)>o 

and clearly the projection of X^{t) on the orthogonal directions with respect to the 
magnetic field oscillates around the magnetic lines. The oscillation frequencies can be 
computed explicitely in this particular case. Observe that the components {Xi,X2) 
satisfy 

xf + ^{{e^) + {9f)xf + ^{{9^} - {effx, = 0. 

The roots of the characteristic polynomial are purely imaginary 

±^y(v^±(^)) 
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d ^ 






~dt ' 




dt 



and therefore the oscillation frequencies in the plane (xi, X2) are ^(a/ {0"^) ± (0)). The 
plasma remains confined along the magnetic lines. Generally we establish the following 
result 

Theorem 2.2 Assume that A G C^(M) is nonin creasing, nonnegative and vanishes on 
[L, +00 [, for some L > 0. Let the initial condition f^^ satisfy 

f p) < X{x{x, p - mUcdiO) /2 ^x) + e0(O, x)) 

where 



1 



2 



2 I I 12 

"^'^c///i2\ /n\2\\ X\ 
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lim |e0(t,x)- rsup{e9,0(.,y)}ds + ^((e2^-(^^)')^l = +oo (16) 



If the electric potential G C^(]R+ x M^) satisfies 

nuu.'"^ |±^|2 

_ , _ ,'e(9t(/)(s,y)| ds + ^ 

uniformly with respect to t E ffi+,X3 G M, t/ien there is a constant R > such that for 
any t G M+, s G M the solution of the problem f ll3p . (1141) verifies 

supp /°(t,s,-,-) C {(x,p) : |^x|<i?}. 

We also prove a strong convergence result 

Theorem 2.3 .4ssiime t/iai E G Li^„^(M+; VT^.ooj^^s^^^ ^ L\^^{^+- L°° {W)) , the 
initial condition /™ has compact support and belongs to iy^'°°(]R^x]R'^). Let f^{t, s, x,p) 
be the solution of f|T3|) . f|T4|) and (/'^)£ i/ie solutions of the problems (I5j), ([2j). Then for 
any interval [0,1] C M+ i/iere zs a constant C{I) such that 

||r(t,-,-)-/°(t,t/£,-,-)||L2(M3xM3) <C(/)£, tG [0,/], £>0. 



3 Average operator 

The previous considerations clearly show that the limit of the Vlasov equation with 
fast oscillating magnetic field deals with multi-scale techniques and homogenization 
arguments. For the rigorous derivation of the limit model (1T2|) we appeal to a slightly 
different method, based on Hilbert expansion at the density level 

nt,x,p) = f{t,t/e,x,p)+ef\t,t/e,x,p) + ... (17) 
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In this section we assume that E = —'Vx4' is a given electric field. Plugging this ansatz 
into (E]) leads to 

dtf + -dsf + e (dtf + -dsf] + ... + ^ ■ (V.r + eVj' + ...) 
e \ e J m 

+ (eE{t,x) + '^e\t/e)^x + u,e{t/E)^p) ■ {y,f + ^^pf + ■■■) = ^ (18) 
and we obtain formally 



Usf ' "^^^-'-^^ 

at the lowest order e'^ and 



dsf + ^e'{s)^x-v,f = o (19) 



dtf + ^ ■ V.r + {eE{t, x) + uj,e{s) ^p) ■ V,f + dsf + -^9 '{s) ^x ■ V,f = (20) 
at the next order e^. The following operator will play a crucial role in our analysis 

ru = div^s,,){u{i,^e'{s)^x)} (21) 

with domain 

D{r) = {ue LliRs-, L\RlxRl)) : div(,,,) [u (l,!^9'{s) ^x)] e 4(R,; L2(R3 xM^))} 

with L'^iRg] X) the space of square integrable T periodic functions m : M — > (X, || • ||x), 
endowed with the norm 

T \ 1/2 

\Hs)\\l ds 



The notation || ■ || stands for the standard norm of L'^(Rs', L'^O^l x I^p)) 

= I / / / \u{s, x,p)\'^ dpdxds \ 

We denote by {S,X,P) = {S, X, P){t; s,x,p) the characteristics of the first order 
differential operator ds + '{s) -^x ■ Vp 

^ = 1 "=0. ^^ = !!^»'(SM)-A-(S(r)) (22) 
or ar ar 2 

with the conditions 

S{0; s,x,p) = s, X{0] s,x,p) = X, P{0] s,x,p) = p. 



It is easily seen that 

S{t; s,x,p) = s + T, X{t; s,x,p) = X, P{t; s, x,p) = p -\ ^{6{s + t) - 6{s)) ^x. 

(23) 

Notice that {x,p — ^^^6{s) -^x} is a complete family of functional independent prime 
integrals of ( 122|) . We introduce the average operator along the characteristic flow (l23l) 
cf. [6] 

1 '■^ 



{u){s,x,p) = 7^ j u{S{t;s,x,p),X{t;s,x,p),P{t;s,x,p)) dr (24) 
= - u[s + r,x,p+'^{e{s + T)-e{s))^x^ dr 

for any function u E L|(]R^; ^^(M^ x M^)). 

Proposition 3.1 The average operator is linear continuous. It coincides with the 
orthogonal projection on the kernel ofT i.e., 

(u) G ker T : / / {u — {u))(f dpdxds = 0, V G kerT. 
Jo JrsJ^s 

Proof. For any function u G L'^(Rs', L'^i^l x I^p)) we have 

\{u)Ws,x,p) <- y \u\^ (^s + T,x,p+'^{e{s + T) -9{s))^x'^ dr 
implying that 

(u) {s) dpdx < — 1 1 1 \u\'^ (s + T,x,p -\ —{6{s + t) — 6{s))'^x^ dpdxdr 

T Jo Jrs Jk3 V 2 / 

= — / \u\'^{s + T,x,p) dpdxdr 

T Jo 

^ • (25) 




T 

Therefore we have 



II {u) II < ||n||, V n G L^(M,; L^(M^. x Mp) 

saying that (■) G £(L|(M,; L2(M3 x RS)), L|(R,; L2(R3 x R^))). It is well known that 
the kernel of T is given by the functions invariant along the characteristics fl23|) 



kerT = I'U G L|(M,;L2(M^ x R^)) : 3 w such that u{s,x,p) = v (x,p-^^^e{s) ^x^] ■ 
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Clearly (u) depends only on x and p — ^^9{s) -^x, cf. flMj) . and thus {u) belongs to 
kerT. Pick a function ip G kerT i.e., 



3 ip : ip{s,x,p) = ip (^x,p ^9{s) -^x^ 



and let us compute / = J^Jj^sJ^i {u — {u))(p dpda;ds. Using the change of coordinates 



q = p 7r~^{^) '^^j foi' fixed (s, x), one gets 



l-T r r 

/ TYIUJ 

{u — {u)) [s,x,q -\ 7r-0{s) ^x ) q) dgdxds 

3 V 2 




'0 Jm:-^ 




^/j{x,q) 



J u(^s,x,q-\ ^9{s) '^x^ ds — T {u)^ dqdx 







and therefore {u) = Proj^err"" ^^r any u G L|l(Ms; -^^^(IR^ x Mp)). In particular (u) = u 
for any u G ker T. □ 

We investigate now the solvability of the equation Tu = v. We have a simple character- 
ization in terms of the kernel of the average operator. Notice that if f = Tu G Range T 
we have for any ip G ker T 

{v — 0)ip dpdxds = / / / Tu (f dpdxds = — u Tip dpdxds = 




'0 JR-^JR-i JO JRiJR-i JO JR-''JRJ 

saying by Proposition 13.11 that {v) = 0. Moreover we have 

Proposition 3.2 The restriction of T to ker (■) is one to one map onto ker(-). Its 
inverse belongs to £(ker (■) ,ker (•)) and we have the Poincare inequality 

\\u\\ < T||Tm|| for any u G -D(T) fl ker (■) . 

Proof. We already know that Range T C ker (■). Assume now that u G D(T) flker (■) 
such that Tu = 0. Since (■) = Proji^^^j.^- we have u = {u) = saying that T|kcr(-) is 
injective. Consider now v G ker (•) and let us prove that there is m G ker (■) fl D(T) 
such that Tu = v. For any a > there is a unique Ua G D(T) such that 

a Ua + TUa = V. (26) 

Indeed it is easily seen that the solutions {ua)a>o are given by 

Ua{s,x,p)= / e°''^v Is + T,x,p —^9{s)^x-\ —^9{s + t)^x] dr. 

JR^ V 2 2 / 

11 



Applying the average operator to fl26l) yields (uq) = for any a > 0. We introduce 
the function 

1/ (r; s, x, pj = y f^s + r, x,p —{)\^s) a; H ^c^l^ + rj xj dr. 

Notice that for any fixed (s,x,p) the function r — )■ V(r; s,x,p) is T periodic, because 
{v) = and thus ||V^(t)|| < T\\v\\ for any r G M. Integrating by parts we obtain 

Ua{s,x,p) = — / e°''^drV dr = / ae°''^V{T; s,x,p) ds 
Jr- Jr- 

implying that 



\\ua\\ < / ae°"||y(r)||dr < T||t;||. 
Jr^ 

Extracting a sequence {an)n such that lim„_>+oo = 0, lim„_>+ooMa„ = u weakly in 
L|(R,; L2(M3 X RS)) we deduce easily that 

ueD{r), ru = v, {u) = 0, ||M||<T||f|| 

saying that (T|kcr(-)) ^ is bounded linear operator and || (T|kcr(-)) ^ ||£(ker(-),ker(-)) < T. 
□ 

Remark 3.1 Notice that T^^ leaves invariant the set of zero average functions, with 
compact support. Indeed, if v & ker (■) has compact support, let us say supp v C 
{(s,x,p) : |x| < i?, \p\ < R} for some R > 0, it is easily seen that for any a > 0, the 
function in (1261) has compact support ( uniformly with respect to a) 

supp Uq, C {(s, a;,p) : |x| < i?, |p| < (1 + m|a;c| ||6'||Loo)i?} 

and therefore the weak limit u = limcx^o 'Wo has compact support. 

Notice that we have the orthogonal decomposition of L^(Ms; L^(]R^. xM^)) into invariant 
functions along the characteristics (12^ and zero average functions u = (u) + (u — (u)), 
since by Proposition 13.11 we have 

(u — {u)) {u) dpdxds = 0. 




'0 



We end this section with the following Sobolev inequality 
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Proposition 3.3 There is a constant C = C{T) such that for any function u G -D(T) 
we have 

lklU^(IR.;L2(iR3xl;3)) < C(r)(||M|| + ||Tm||). 

In particular for any function u G D{T) fl ker (■) we have 

< C{T){l + T)\\ru\\. 

Proof. Without loss of generality we can assume that the function u is smooth (the 
general case follows by standard density arguments). For any s G M and t G [s — T, s] 
we have 

\(s,x,p- ^(^(t) - e{s)) ^x] I = (Tn) fs, x,p - ^(^(t) - e{s)) 



ds r V ' 2 ' " yj ^ ' V ' 2 

After integration one gets 

u (^s,x,p-^m - 9{s)) ^x) = u{t,x,p) + j\u (^r,x,p-^m - 9{t)) ^x) 
implying that 

W'^i^y ■) ■)IIl2(r3x]r3) < '2\\u(t, ■, ■)IIl2(r3x]r3) + 2T\\Tu\\'^ . 

Averaging with respect to t over the period [s — T, s] yields 

2 2 2 2 

II'"IIl^(IRs;L2{r3xR3)) ^ +2r||T'u|| 



and the first statement follows with C(T) = ma,x{y^2/T, \/2T}. Moreover, if u G 



D{T) n ker (■) we know by Proposition 13.21 that ||n|| < T||Tm|| and therefore 
IMIl-(r.;L2(m3xr3)) < C{T){T\\ru\\ + lir^ill) = C{T){1 + T)||ru||. 



□ 



4 Commutation properties of the average with re- 
spect to derivations 

We have seen that T^^ restricted to zero average functions is linear and continuous. In 
view of further regularity that we need for the asymptotic analysis of (|5]) we investigate 
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now the action of T^^ and (■) over subspaces of smooth functions. We formulate our 
statements in the general framework of a characteristic flow associated to smooth, 
divergence free fields. Let b : M™ — )■ be a field satisfying 

6 G W^i'^f (M'"), div^6 = 

and the growth condition 

3C>0 : |%)|<C(l + |y|), yeR"\ 
Under the above hypotheses the characteristic flow Y = Y{s; y) is well defined 

— = biYis;y)), is,y)ERxR^ 
YiO;y) = y, y G M'", 

and has the regularity Y G x M™). Since the field is divergence free, we deduce 

by Liouville's theorem that for any s G M, the map y — > Y{s; y) is measure preserving. 
Notice that we don't make any periodicity assumption on the flow Y. As usual the 
notation b ■ Vy stands for the first order differential operator with domain 

D(6 ■ V,) = {ue L'iW^) : dwyiuiy)biy)) G ^^(M™)} 

which maps any u G D(6 ■ V^,) to the function dwy{u{y)b{y)). The kernel of this 
operator is given by functions which are constant along the flow Y 

ker(6 • Vy) = {ue L'^{W) : u{Y{s; y)) = u{y), s G M, a.e. y G M™}. 

It is easily seen that for any function u G L^(]R™) and any T > 



^l\{Y{s;-))ds 



< I|m||l2(ir™) 

L2(IR™) 



and therefore we can expect compactness properties for the family of averages along 
the flow Y i.e., {j; JqU{Y{s;-)) ds,T > 0}. Indeed, the mean ergodic theorem, or 
von Neumann's ergodic theorem (see [18], pp. 57) allows us to construct the average 
operator along the flow Y. 

Proposition 4.1 For any function u G the averages |t u{Y{s] ■)) ds con- 

verge strongly in when T — )■ +oo, towards some function denoted (u) G 
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L^(M'"). The map (•) : L^(]R™') — )■ L^(R™) zs linear, continuous and coincides with 
the orthogonal projection on ker(6 ■ V^) 

{u) e ker(6 -Vy) : [ {u{y) - {u) {y)My) dy = 0,W^e ker(6 ■ V^). 

It is easily seen that Range (b-Vy) C ker (■). Actually it is shown in [6] that Range (6 ■ Vy) = 
ker(-). We are searching now for derivations commuting with the average operator. 
We prove that the average operator is commuting with any derivation c ■ Vy associated 
to a field c in involution with b. We recall here the following basic results concern- 
ing derivation operators along fields in M"*. For any C, = {C,i{y), ■■■,C,m{y)), where 
y G M™, we denote by the operator ^ ■ \/y. A direct computation shows that for any 
smooth fields ^,ri, the commutator between L^,L^ is still a first order operator, given 
by [L^, Lr,] := L^Lf^ — L^L^^ = L^, where x is the Poisson bracket of ^ and rj 

X=[^,V], [^,v]i = {^-^y)Vi-{v-^y)^i = L^{r]i)-L^{^^), i e {1, ...,m}. 

It is well known (see pj, pp. 93) that L^,Lfj commute (or equivalently the Poisson 
bracket [^,r7] vanishes) iff the flows corresponding to ^,r7, let say Zi, Z2, commute 

Zi{si;Z2{s2;y)) = Z2{s2;Z,{s,;y)), s,,S2eR, y G M"^. 

Consider a smooth fleld c in involution with b and having bounded divergence 

c e VTiof (M™), divj^c G L°°(M'"), [c, 6] = 

and let us denote by Z the flow associated to c (we assume that Z is well defined for 
any {s,y) G M x M"*). We claim that the following commutation property holds true. 

Proposition 4.2 Assume that c is a smooth field in involution with b, with bounded 
divergence and well defined flow. Then the average operator commutes with the trans- 
lations along the flow of c 

{uo Z{h;-)) = {u) o Z{h;-), ueL\W"), /i G M. 
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Proof. The commutation property of the flows Y, Z and Proposition 14.11 allow us to 
write the strong convergences in L^(R'^) 



{uoZ{h;-)) 



(27) 



1 r 

hm - / uo Z{h;Y{s;-)) ds 

T-^+oo T Jq 

1 /-^ 

hm - / uoY{s;Z{h;-)) ds 

T^+OO T Jq 

(m)oZ(/i;-). (29) 



Notice that the third equality in the above computations follows by changing the 
variable along the flow Z and by using the boundedness of div^^c. ^ 

We denote by c ■ Vy the operator with the domain 

D(c- V,) = {uE L^W^) : dwy{u{y)ciy)) E L\R"')} 
which maps any function u G D(c ■ Vy) to the function divy{uc) — u div^c. 
Proposition 4.3 Under the hypotheses of Proposition assume that u G D(c- Vj^). 



Then {u) G D(c ■ Vy) and c ■ Vy {u) = (c ■ Vyu) . 

Proof. For any /i G M* we have 

{u) o Z{h; ■) — {u) {u o Z(h; •)) — {u) / uo Z{h; ■) — u 



h h \ h 

Since u G D(c ■ Vj^) we have the strong convergence \imh-^o{u o Z{h; ■) —u)/h = c - VyU 
in L^(]R™) and by the continuity of the average operator, we deduce that 

(u) o Zih: ■) — (u) , „ , 
hm ^ ^ = c ■ Vyu) 

strongly in L^(M'"), saying that {u) G D(c ■ Vy) and c ■ Vy {u) = (c ■ Vyu). □ 
Similarly we can prove 

Proposition 4.4 Assume that u G iy^'^([0, T]; L^(]R™')) for some p G (1,+cx)). Then 
the application {t, y) — )■ {u{t, ■)) {y) belongs to iy^'P([0, T]; L^(M™)) and we have dt (u) = 
{dtu). 
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We are ready now to study how the regularity propagates under the action of {b-Vy) ^. 
We make the following assumption (Poincare inequality) 

3 Cp > such that ||'u||L2(igm) < Cp\\b ■ Vyu\\L2(^m^, V m G D(6 • Vy) fl ker (■) (30) 



meaning that the range oib-Vy is closed i.e., Range (6- Vj^) = Range (6 ■ Vy) = ker (■) 
and b ■ Vy restricted to ker (■) is one to one map onto ker (■) with bounded inverse. 
Notice that the above hypothesis is satisfied by the operator in fl^T]) . cf. Proposition 



Proposition 4.5 Under the hypotheses of Proposition o.nd fl30|) assume that v G 
ker (•) n D(c ■ Vy). Let us denote by u the unique zero average solution ofb ■ VyU = v. 
Then u G D(c ■ Vy) and \\c ■ V yu\\i2(^^ra'^ < Cp\\c ■ V^f ||i2(igm). 

Proof. For any /i G M* we have Uh '■= uoZ{h; ■) G D(6- Vj^) and b-VyUh = {b-VyU)h = 
Vh. Therefore we deduce that b ■ Vyiu^ — u) = v^ — v. Since the average is commuting 
with the translations along the flow of c we have 

{uh -u) = (uh) - (u) = {u)^ - (m) = - = 

and we can apply the Poincare inequality (1501) 

\\Uh — < Cp\\Vh — f ||L2(|5m), h G M*. 

Therefore sup^^jg* \\uh — w||l2(M"i)/|/''| < Cp\\c • V yv\\L2^^m^ saying that u G D(c ■ Vy) 
and ||c ■ VyM||L2(]gm) < Cp\\c ■ Vj,f ||z,2(]Rm). Notice also that we have 

b ■ Vy{c ■ VyU) = C ■ Vy(6 " V yU) = C " V yV 

with (c ■ Vyu) = c - Vy (u) = saying that (6 ■ Vy)^^{c ■ Vyv) = c ■ VyU. ^ 



4.1 Regularity propagation under fast oscillating magnetic 
fields 

We apply the previous general results to the operator fl?I]) acting on the phase space 
(s,x,p) G R^. We indicate a complete family of fields in involution with respect to 
ds + ^^^6'{s) -""X ■ Vp. The reader can convince himself by direct computations. 
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Proposition 4.6 The following fields are in involution with respect to '(s) -"-a;- 

• '^{s,x,p) = ^(^{s)dp2, (? ■ '^{s,x,p) = ~^ — ^^(■5)<9pi, ■ '^(s,x,p) = 9^3 

C ■ ^ {s,X,p) 9piy C • ^ (^s,X,p) 9p2y C • ^ [s,X,p) • 

Proposition 4.7 Assume that v G ker (■) such that VxV^VpV G (L|^(Rs; x 
Mp)))'^. Let us denote by u the unique zero average solution of Tu = v. Then 
VxU,VpU G (L|(M,;L2(M3 X R^))^ and 

\\dxiu\\ < T{\\dccM\ + "^|wc|||6'||l°°(r3)||<9p2^II} 

\\dx2u\\ < T{\\dcc2v\\ + m\uJc\\\0\\L^(^M3)\\dpM\} 
\\dxM < T\\dxA, \\dpM\ < T\\dpM\, ^ e {1,2,3}. 

Proof. By the hypotheses we know that v G nf^^D(c* ■ V(s,x,p)) and therefore Propo- 
sition 14.51 implies 

u G ntiD(c^ ■ V(,,^,p)), \\d ■ V(,,^,p)u|| < T\\d ■ V(,,^,p)^;||, i G {1, 6}. 
In particular we have V^u, VpU G (L|(]R,; ^^(IRS x Rl))f 

\\dxA\<T\\dxAl \\dpM\\<T\\dpM\, ^€{1,2,3} 

II II 1 V7 , ^^Cq r. II 

\\Ox^u\\ = \\c ■ V(s,x,p)U + -^0 dp^uW 

Tfl OJ 

< \\C^ ■ V(s,^,p)M|| H Y^\\0\\L°-{R3)\\dp2u\\ 

< T\\c' ■ V^s,x,p)V\\ + T^^\\e\\L^^^s)\\dp2v\\ 

< T{\\d,M\+M^c\\\0\\Lo^iUS)\\dp,v\\} 

and similarly ||<9^.2m|| < T{||9^2i;|| + m\uJc\\\9\\L--(R^)\\dpM\}- □ 
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5 Limit model 



We are ready now to investigate the limit model of ([5]) as e \ 0. We appeal to the 
method introduced in [6] for general transport problems (see also [4J) which combines 
Hilbert expansion and average properties. A rigorous convergence result is obtained 
using the notion of two-scale convergence. We analyze the properties of the limit model, 
in particular we establish the energy conservation and justify the confinement around 
the magnetic lines. The terms in the Hilbert expansion (IT7|) satisfy 

rf = (31) 

dtf + - ■ V./° + {eE{t, x) + u,9{s) ^p) ■ Vpf + Tf = 0. (32) 

m 

The equation (13T1) says that at any time t G R+ the function {s,x,p) — )■ f^(t, s,x,p) 
belongs to ker T and therefore 

nt,s,x,p) = g' (t,x,p-^^(.) ^x). 

The time evolution equation for /° is obtained from fl32|) after eliminating Thanks 
to Proposition 13.21 we have for any t G M+ 

dtf + - ■ V,./° + {eE{t, x) + oo,e{s) ^p) ■ V^f G Range T = ker (■) . 
m 

Therefore ( 132|) is equivalent to 

(dtf + - ■ V,/° + {eE{t, x) + u,9{s) ^p) ■ V^f) = 0, te M+. (33) 
\ m I 

It is easily seen that {dtf)=dt{f)=dtf. It remains to compute the averages of 
the derivatives with respect to x and p. For simplifying the computations we assume 
that /° is smooth but it can be shown that the limit model which we will obtain still 
holds true in the distribution sense. 

Lemma 5.1 Assume that f{s,x,p) = g {x, q = p — ^^9{s) -^x) is smooth. Then we 
have 

= - y - (0)) ^x) . V./ + (y - (0)) + ^{29 {9) - (9^) - 9^) ^^x) 
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and 

= [eE + {9) + ^{{9') - 9 {9)) ^^x^ ■ V,/. 
Proof. We have, with the notation -""Vg = {dq^, —dq-^,0) 

TYILl) 

VJ = V.g + ^9{s) ^Vqg 

and therefore 

■ V./) = (P) . ^ + f (pfl) . -V,,, 

\m I ml 

since the derivatives of g are constant along the characteristic flow fl22l) . By the defi- 
nition of the average operator we have 

(p) = - y |p-_^^(s)^a;+^^(r)^a;| dr 
, "^^^ ia\ ^ 

where q = p ^9{s) '^x. Similarly we obtain 

{p9) = ^£ {p-^9{s)^x + '^9{r)^x}9{r)dr 

= q{0) + ^ {e') 

and finally 

The second statement follows easily observing that Vp/ = Vqg and therefore 

{{eE{x)+u,9^p)-Vpf) = eE{x)-Vqg + Uc{9^p)-Vqg 

= eEix) ■ Wqg + [uj, {9) + ^ {9') ^^x^ ■ W^g. 

□ 

Combining the previous computations yields the transport equation 
dtg' + (I + y {0) ^x) . V./ + (ei? + y {0) + ^ {e') ^^x^ ■ V,^?° = (34) 
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which is exactly f|T2|) . Performing the change of unknown f^[t, s,x,p) = g^{t,x,p — 



^e{s) ^x) leads to 



dtf + {^-^{e{s) - {e))^x) - VJ' (35) 

+ (eE{t^ x) + ^{e{s) + {6)) + "^{{e') - 0\s)) ^^x) . V,r = 0. 

Based on the concept of two-scale convergence, introduced in [17\ and developed in [1] 
we prove a two-scale convergence result of the solutions in (jS]) towards (^S^, supple- 
mented by an appropriate initial condition. 

Definition 5.1 Let (/"(t, be a sequence in L^([0, /]; L^(M^ x M^)). We say 

that {f^)n two-scale converges towards some function /° G L^([0, /]; L|t(Rs; L^(Ri^. x 
Mp))) iff we have 

lim / / / f^{t,x,p)ip{t,t/e,x,p) dpdxdt = [ ^ [ [ [ {f°ip){t, s,x,p) dpdxdsdt 

n^+oo Jq J^3jf,3 Jq I Jo Ju^JuS 

for any test function (p G ^^([0, /]; (R,; ^^(RS x R^))) (here (M,;X) stands for 
the set of continuous T periodic functions with values in the normed linear space X). 

Adapting the arguments in we obtain 

Proposition 5.1 Let {f'^{t,x,p))e>o be a bounded family in L^([0, /]; L^(]R^ x M^)). 
Then there is a sequence En \ and a function /° G L^([0, /]; L|^(Rs; L^(]R^ x Mp))) 
such that {f'^")n two-scale converges towards /°. 

Proof, (of Theorem 12. ip We use the weak formulation of ([5]) with the test function 
ri{t)(p{t/e,x,p) where r] G Cl{R+) and ^ G Cl{R/TZ xW x R^). We obtain 





V / / f{t,x,p) 
4. Jr^Jrs 



1 p 

-dsf(t/e,x,p) H Vx'~p(t/e,x,p) 

e m 



^(0) / / /'"(a:;,p)</'(0, dpdx — / ri'{t) j j f^{t,x,p)(p{t/e,x,p) dpdxdt 

dpdxdt 

- V / f%t,x,p) \eE +^^e'{t/e) ^x + Uc9{t/e) ^p] ■ Vp<p{t/e,x,p) dpdxdt 
= 0. 

Multiplying the previous formulation by e, one gets by two-scale convergence (after 
extraction eventually) 

V{t)-k I I I f°{t,s,x,p)Tf{s,x,p)dpdxdsdt = 
^ Jo Jr^Jr^ 
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saying that f^{t, ■,-,■) G ker T for any t G ^4.. We use now the same weak formulation, 
but with (p G Cl{R/TZ x x M^) fi kerT. For example take (p{s,x,p) = ilj{x,p - 
^e{s) ^x) with G C1(R3 X M:^), let us say supp C : max{|x|, |g|} < R}. 

In this case supp (y9 C {{s,x,p) : max{|x|, < (1 + mlwd ||6'||loo(]r3)/2)_R}. Since 
ip{s,x,p) G ker T we have 

1 TTlLfJ 

-ds(f{t/€, X, p) + —^6' '(t/e) ^x ■ Vpip{t/e, x,p) = 
and the weak formulation becomes 



—77(0) / / f^{x,p)(p{0^ x,p) dpdx — / ri'{t) / / f'^{t^x^p)(p{t/e,x,p) dpdxdt 





p_ 

m 



Vit) / f -■ VMt/^, X, P) + {eE + Uce{t/e) ^p) ■ VMt/e, x, p) 



0. 



dpdxdt 
(36) 



As before, by two-scale convergence one gets 



lim / f]'{t) I I f'^{t,x,p)(p{t/e,x,p) dpdxdt 

Jr+ jR^ijM^i 

r]\t)^ I I I f{t,s,x,p)<p{s,x,p)dpdxdsdt. 

J Jo jRsJ^a 



Similarly one gets 




\im I r] I If—- VMt/e,x,p) + (eE + uj,e{t/e) ^p) ■ Vp^{t/e,x,p) 




fit, s, x,p) -■ V,<^ + {eE + UcOis) ^p) ■ 



1 



dpdxdt 
dpdxdsdt. 



Therefore fl36|) implies 



—77(0) / / f^{x,p)(p{0, x,p) dpdx 




s, X, p)^{s, X, p) dpdxdsdt 



T 




Jr^J«.3 



T Jo Jr^Jrs 

fit, s, x,p) — ■ V^v? + ieE + UcOis) ^p) ■ Vpf dpdxdsdt. (37) 
m 



We transform the term involving the initial condition. Since (f> G ker T we have for any 
s G M 

(^(0,x,p) = (^s,x,p ^^(0) 3;+-^^(s) xj 
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implying that 



tI [ fM^,x,p)dpdx = [ [ [ fix, p)y^(s,x,p-^{e{o)-e{s))^x) dpdxds 

f [x,p- ^{9{s) -9(0)) ^xj ipis,x,p) dpdxds. 




(38) 



We transform now the right hand side of fl37j) thanks to Lemma 15.11 Indeed, since 
f{t, ■,-,■) G kerT, for any t G M_|_, we have by the variational formulation of the 
average operator in Proposition 13.11 

f{t) \— ■ V^if + {eE + Uc6is) ^p) ■ Vp J dpdxds 




'0 Jm?. 



m 




'0 



r-T 




+ ( + |(^(s) + {6)) ^p + "^{{6^) - e\s)) ^^x ) ■ Vp^ 



f {—■ V^v^ + ieE + uOcOis) ^p) ■ V„V9 ) dpdxds 
i \m ^ 'I 

dpdxds. (39) 

Combining now (l37j) . (I38i) . (139|) implies that /° is a weak solution of the transport 
problem (IT^ . (fill) . Equivalently, the function such that /°(t,s,x,p) = (?°(t,x,p — 
mcJc/2 6'(s) -""x) is a weak solution of the transport equation supplemented by the 
initial condition 

^7°(0, X, g) = r (x, g + ^^(0) ^x) . (40) 

□ 

Notice that the model (1341) . (1401) is posed in a six dimensional phase space whereas 
the model f|T3l) . f|T^ acts on a seven dimensional phase space. Thus, at least for the 
numerical point of view it is preferable to appeal to f lM|) . fHOj) . 

It is interesting to observe that, as e \ 0, the kinetic energy of /°(t, t/e, ■, ■) ~ 
/^(t, - , ■) can be decomposed into kinetic energy and elastic energy associated to the 
density 
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V^\t) = IJJyi^ ^x" + ^((^^>-(^)^)^j,°(t,x,g)dgdx. 

Moreover, when the electric potential solves the Poisson equation corresponding to the 
concentration 

/ /°(t,t/£,x,p)dp= \ g''(t,x,p-'^e{t/e)^x) dp= I g'{t,x,q)dq 
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2 I I 12 

muj^ , ,„o\ 2x x' 




the total energy is preserved. In order to simplify our computations we work with 
smooth solutions. 

Proposition 5.2 Assume that is nonnegative, e 1^^'°°(R^ x B?) such that 

(1 + \^x\^ + \p\^)f''{x,p) dpdx < +00. 
1 

iJIfEe i^ioc(K+; Vr^'°°(M'^)) is a given electric field then 

[ [ ^-^f{t,t/e,x,p) dpdx ^W^{t) weakly ★ in L^^{R+) as e\0. (41) 

a) If f^'^ has compact support and E e Lj'Qj,(M_|_; W^'^ i^)) , then there is a continuous 
nondecreasing function R depending on t, such that 

supp^°(t, •,•) C : \x\<R{t), \q\<R{t)], t G M+. 

Hi) If E belongs to e Lj'-pp(R+; VF^'°°(]R^)) and solves the Poisson equation 

E^-W^(P, -A^(P{t)^-( I g\t,x,q)dq-no{x)] 

such that E(0, •) G {p{^^)f, then 

dt \ 2 7^3' ' J 

Proof, i) It is easily seen that 

/ / f^{t,t/e,x,p) dpdx = 9^ (t,x,p — '^^^^^^Oit/e) '^xj dpdx 

= / / 9°{t,x,q) dqdx 
Jr3 Jr3 

= 11 g\Q,x,q)dqdx 

Jr» JR3 



(42) 



/ / r(x,g + ^^(0)^x) dqdx 




f^^{x,p) dpdx 
i 

saying that f^{t,t/e, ■) belongs to L^{Rl x R^) uniformly in i e R+, e > 0. We 
consider the function 



1 r 1 



T 



2m 



<l+^0{s) ^X 



ds 



1 

2m 



4 

moo, 
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Observe that 

+ T <"> ('^ + f (9) + =^ <«^> ^^x) .V,x = e£. + f (9) 

(43) 

and thus multiplying fjM|) by x we deduce 



/ [g\{x,q)dqdx = [ f g\t,x,q)eE ■ + {6) ^x) dqdx (44) 




< I I g'\eE\\ — q+'^{9) ^x 



1 



2m 



2 



2m 



dgdx 






< ||eE(t)|U.(,3) / / , V dgdx + Mi^^J^ / I rdpdx. 

By Gronwall's lemma it follows that /jga/j^g x, g) dgdx G i^i^c(^ 

^°(t,x,g)xdgdx < (! [ ^°(0,x,g)xdgdx+ "'^""^fJ^"°°^^^» / / r dpdx 

X exp ||eE(s)||Loo(]R3) ds 

For any function t] G Lj'Q^(R_|_) we can write 

lim / 7]{t) [ [ f{t,t/E,x,p)^-^dpdxdt 
''\^Jr+ JiR3JiR3 2m 

= hm / vit) [ [ g' (t,x,p- ^^eit/e) ^x) ^-f- dpdxdt 



^N'O Jr^ iiR3j]R3 2m 



lim / ri{t) I I g^{t,x,q)^ q + ^^0{t/e)^x 



2 



dgdxdt 




Vit) / / 9°{t,x,q)x{x,q) dqdxdt 




saying that 

^-^f\t,t/e,x,p)dpdx^W^{t) weakly * in Li^,(R+) as e \ 0. 
s2m 

ii) Assume that supp /™ C {{x,p) : |x| < -R™, |p| < i?™} with i?™ > 0. Therefore 
supp(7°(0,-,-) C {(a;,g) : |x| < \q\ < {l + m\uj,\/2 \\e\\L^^^^)) R'^} 



C {(x,g) : v6K^<^^^",k3| <^i?"^} (45) 
for some u > 0. We consider a function ,^ satisfying 

^eC\R), 0<e<l, suppe = K+, r>0. 
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Applying the weak formulation of flMl) with the test function (t, x, q) — )■ ^{^/xi^^ q) ~ 
a{t)), with a E C\R+) yields by 



g'^ii, X, q)^i\/ xix, q) - a{t)) dqdx - / / g^{Q, x, q)i{\/x{.x, q) - a(0)) dqdx 

eE{T, x) f q (jJc 




g\T,x,q)C {-a'{r) + 



< f f /«V...,)r|-^'(r)+ ''^'^''''"il-'''' |d,d.dr. 



f ^ + ^ (0) I dgdxdr 
Vm 2 / J 



We take a(0) = uR"" and a'(r) = ""^^"J^"""'' , r G M+. Since g%0,x,q)^{^/x{x,q) 
q;(0)) = one gets 




g^ii, X, q)^iVxix,q) - a{t)) dqdx < 



implying that supp g^{t,-,-) C {{x,q) : a/ x(a;, g) < a(t)}. Similarly, applying the weak 
formulation of (134|) with the test function (t,x,q) — > ^(|x3 — tga/ml — /3 G C^(R+) 
we obtain 




g^{t,x,q)^{\x3-tq3/m\-l3{t))dqdx- / / (?"(0, x, g)e(|x3| - /3(0)) dgdx 
= f f /" |-/3'(r) - ^^^^^^^^sgn(x3 - tgg/m) j dgdxdr 

Jo JR3yR3 I rn J 



Taking /3(0) = z/i?'" and /3'(r) = r G M+ we deduce that 

supp ^(^(t, ■) C : |x3 - tga/m] < 

We have proved that 



supp5(°(t) C {(x,g) : < J^^" + ^==II^IUM[o,t];J^-(iK3))} 

V2m 

n {{x,q) : |x3|<z^i?- + tM + '"'^'"^'"°'*^^"° 



m m 
C {(x,g) : |a;| < R{t), \q\ < R{t)} 

for some continuous nondecreasing function R{t). Notice also that for any t G M+ 

supp •,-,■) C : \x\<R{t), IpI < (l + m|u;,|/2 ||^|Uoo(K3))i?(t)}. 
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iii) Assume now that the electric potential solves the Poisson equation fj42|) . By stan- 
dard computations involving the continuity equation 

dJ g'dq + divj (^ + ^{9) ^x)gUq = 



one gets 



d Eq 



[ \E\^dx = - [ eE{t,x)- [ (- + '^{6) ^x)g'' dqdx. (46) 



Combining (04]), (gH]) yields 



d 
dt 



[ [ g\t,x,q)x{x,q)dqdx + ^ f \E{t,x)\^ dx] = 0. 

JR3JM3 Z J ^3 ) 



□ 

In the sequel we inquire about the confinement properties of the limit model fl5^ . f HU]) 
(resp. ffT^ . dH])). We exploit here the invariants of fl51|) . Let us consider for the 
moment that the potential is stationary. In this case notice that (IMI) writes under the 
Hamiltonian form 

dtg' + V,H ■ V,^° - V,H ■ V,(7° = 

with the Hamiltonian 



1 



H{x, q) = xix, q) + e0(x) = ^ q + {9) ^x' + '-^{{9^) - (9)')^ + e0(x). 

In particular H is a. stationary solution of or equivalently H is an invariant of the 
characteristic flow of (13^ 



2 "^I^c///i2\ //i\2x|"'"^P 



^ = v,H{x{t), g(t)), ^ = -v.if(x(t), g(t)). 

Under additional hypotheses on the electric potential it is easily seen that the plasma 
remains conflned in a bounded region around the magnetic fleld lines. Indeed, assume 



that the hypothesis in Theorem 12.21 holds true 

hm {e<P{x) + ^{{9^)-{9f)^-^] = +oo (47) 
uniformly with respect to X3 and that at the initial time we have 

< g^{0,x,q) < X{H{x,q)), {x,q)eR^xR^ (48) 
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for some nonnegative profile A G C^(M), vanishing on [L, +oo[. By the maximum 
principle we deduce that 



< g%t, X, q) < X{H{x, q)), {t, x, g) G M+ x x 

which guarantees the compactness of the support of g^{t) along the orthogonal direc- 
tions to the magnetic lines, uniformly in time. Indeed, we have 

supp^°(t)c{(x,g) : Hix, q) < L} C {{x, q) : e<Pix) + "^{{6') - (9)')^-^ < L}. 

By the hypothesis (l47|l there is i? > such that for any x G verifying \-^x\ > R we 
have 

e<^ix) + ^i{e^)-{e)i-^>L 

and finally 

supp ^°(t) C : \^x\ < R}, VtGM+. 

Remark 5.1 // the electric potential depends only on {-'-xl and x^ then {q ■ -""x) is 
another invariant of fl34|) . // the electric potential depends only on [-""xl then qs is an 
invariant of f p4|) as well. 



When the electric potential depends on time, the previous Hamiltonian becomes H{t, x, q) 
x(x, q) + e(j)(t, x) and therefore we obtain 



c /n2\ ±±. 



^ ^9') ^^x ] ■ V,H = edtct). 



\m 2 / \ 2 

(49) 

In this case we need to construct a particular super-solution for 



Proof, (of Theorem 12.21 for time dependent electric potential) We work in the phase 
space (x,g), using the equation and the initial condition fHOj) . By the hypotheses 
we know that 

^°(0, X, q) < X{x{x, q) + e0(O, x)), (x, g) G x 
Let us consider the function h defined for any (t, x, q) G R+ x x by 

h{t, X, g) = A ( x(x, g) + e0(t, x) - [ sup {edt(p{s, y)} ds ] . 

28 







Observe that 

h{0, X, q) = X{x{x, q) + e0(O, x)) > g°{0, x, q), {x, g) G x Ml 

Taking into account it is easily seen, by the monotonicity of A, that 

d,h + + y (0) ^x) . V^h +{eE+^ (6) ^q + ^ {6') ^^x^ ■ V ,h > 0. 

By the maximum principle we deduce that 

g°{t, X, q) < hit, X, g), (t, x, g) G M+ x x 

By the hypothesis (|T6l) there is i? > such that for any (t, x) G M+ x M?, with [-"-xl > i? 
we have 

2 I _L 1 2 

e0(t,x)- / sup{e9,0(.,i/)}d. + ^((^^2)-(^^)^)^>L 
implying that for any t G M+ 

supp5'°(t) C {(x,g) : x(x, g) + e0(t, x) - / sup{e(9t0(s, y)} ds < L} 

Jo weK3 



/'t 2 I _L 1 2 

/ sup{e9,0(., y)] ds + "^{{e') - {efY-^ < L] 
Jo weK3 4 z 



C {(x,g) : e0(t,x) 
C {(x,g) : |^x| < R}. 

Since f'^it, s,x,p) = g^(t,x,p — mUc9{s)/2 -^x) finally one gets 

supp /°(t,s, -,■) C {(x,p) : |^x|<i?}, tGM+,sGM. 

In particular 

supp -,■) C {(x,p) : |^x| < i?}, t G M+, e > 0. 



□ 



6 Asymptotic behaviour 

The aim of this section is to justify rigorously the Hilbert expansion f[T7|) . More 
precisely we intend to prove that /^(t, x,p) = t/e, x,p) + 0{e) strongly in L^(M'^ x 
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M^), uniformly for t in bounded intervals. The idea is to introduce the solution F"^ = 
F^{t, s,x,p) of the transport problem 

dtF' + — ■ V^F' + (eE(t, x) + u^eis) ^p) ■ V„F' + -TF' = (50) 

m e 

with the initial condition 

F%0, s, x,p) = r - ^e{s) ^x + ^^(0) ^x) (51) 

and to observe that F'^(t,t/e,x,p) satisfies (j5]), (EJ, saying that F'^it^t/ e^x^p) = 
f^(t, x,p). We start by estimating the error between F^(t, s, x,p) and /°(t, s, x,p). We 
also prove that for any / G M+ the functions {F^{t) : t G [0,/],e > 0} are uniformly 
compactly supported. 

Proposition 6.1 Assume that E G Lii„,(M+; VT^.ooj^^S))^ g^^ ^ Li^^(M+; L°°(M3)) for 
any I G M+, /™ G iy^'°^(M^ x M^) and supp /™ compact. Then for any I G M+, 
i/iere is a constant Ci{I) such that 

\\F'it)-fit)\\<C^iI)e, tel, E>0. 

Proof. By Proposition 15.21 ii). the solution of the problem f[T^ . f lT^ has compact 
support, uniformly for t in bounded intervals i.e., V / G M+, 3 R{I) such that 

supp/°(t,-,-,-)c{(s,x,p) : |x|<i?(J), \p\<R{I)}, te[0,I]. 

Under our hypotheses, the solution /° has the regularity 

Recall that the model (135|) is equivalent to (133|) . By Proposition 13. 2[ for any t G M+, 
there is a unique function f^(t) of zero average such that 

dtf + £ ■ V,/° + (eE(t, x) + u,e{s) ^p) ■ Vp/° + rf\t) = 0. (52) 

Since /° is smooth and has compact support (uniformly with respect to t in bounded 
intervals) the following set 

{V(t,,,p) (dtf + £ ■ V,/° + (eE(t, x) + u;e^(s) ^p) ■ Vpf) : t G [0, /]} 
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is bounded in L|(Ms; L'^i^l x M^)). By Propositions lOl SiZl we deduce that 

{V(t,.,,)/^(t) : te[0,I]} 

remains bounded in L'^CRg] L'^ {M-l x Mp)). Combining (150|) . (152|) and the constraint 
dSH) yields 

(5* + ^ ■ v. + {eEit, x) + ^p) ■ V,) {F^ -f- ef} + - - e/^} 

= -e {dtf^ + ^ ■ V./^ + (eE(t, x) + u,e{s) ^p) ■ Vp/^l 
I m J 

and after integration with respect to one gets 



d 



^\\F'-f- II < 5 dtf + ^- VJ' + {eE{t) + co^e ^p) ■ VJ' 
at m 



Taking into account that F^, / satisfy the same initial condition we deduce that 
\\F^{t)-f{tyef\t)\\ < s\\fm+e / + ^ ■ V. + {eE{r) + oo^O ^p) ■ V,) 



dr 



and finally 



|i^^(t)-rWII <Ci(/)£, ie/, e>0. 



□ 



Proposition 6.2 Assume that the electric field is smooth E G iy-^'°°(]R^)) 
and that /™ has compact support supp /™ C : |x| < i?™, |p| < T/ien it/iere 

zs a continuous nondecreasing function 6 : — t- R_|_ swc/i i/iai 

suppF^(t,-,-,-)c{(s,x,p) : IpI < 5{t), "^l^-NI^II^- |a;| < 5{t)}, teR+,e>0. 

In particular for any I G IR+ 

U suppF^(t,-,-,-)c{(s,x,p) : \p\<5il), !I^MMl^\^\<s(^i)y 
£>o,te[o,/] 

Proof. Let us consider a function ^ G (M) satisfying 



o<e<i, r>o, suppe 



for some constant C > 0. We denote by the function 

m\uc\ 



h{s,x,p) = p ^6{s) '^x 



+ 



\0\\ TOO |x| 
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Notice that h depends only on the invariants x and q = p — ^^^6{s) -^x and therefore 
Th = 0. By direct computations one gets 



and therefore 



P ^ 



p U[S) X 



q_ 



P -^e{s)^x = eE{t)-f- + ^—-^^^x ^ 



\q\ 

< ||eE(t)|Uco(K3) + 



4 



mu 



\q\ 



c II/3II2 



\X\. 



Similarly 



^ ■ V. + {eE{t, x) + u,d{s) ^p) ■ Vp ^ ^ 
m J I 



\x\ 



p m\ujf\ \Q\l°° X 
m 2 



x\ 



< 



Finally the function h satisfies 

(dt + --V, + ieE{t,x)+Uce{s)^p)-Vp + -r]h < ||e^(t)|Uoc(K3) + ^ 
\ m e J 

Using now the weak formulation of f lSU]) with the test function (t, s, — )■ C,{h{s, x,p) — 
6{t)), with 6 e C\R+) yields 



-h. 



rT r r rT 

F'{t)^{h{s,x,p) - 5{t)) dpdxds - 




'0 jm»Jr3 




JR^Jm? 



F%0)^{h{s,x,p) -5(0)) dpdxds 



< 




F^'{T)C{h{s,x,p) - 5(r)) -5'(r) + ||eE(r)|Uoo + 



'0^0 



h f dj9da;dsdr. 

(53) 



Notice that 



supp F^(0, ■,■,■) C {(s,x,p) : |x|<i?^ b|<(l + m|a;,| \\9\\loo)R''} 

and therefore 6o := sup{h{s,x,p) : (s,x,p) G Ue>oSupp F^(0)} < +oo. In this case we 
have 




^'^(O, s, x,p)^{h{s, x,p) — Sq) dpdxds = 0, e > 0. 



(54) 



We determine the function 6 by solving 



5\r) = \\eEir)\U^+\-^^^^^ re 
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with the initial condition 6{0) = 6q. The right hand side of fl53|) becomes 




F%T)C{h{s,x,p) - 11^11^°° (/^ _ S{t)) dpda;dsdr 



'0^0 



< c^^^^^^J^ [ [ [ [ F%T)^{h{s,x,p) - S{t)) dpdxdsdr. (55) 

2 JoJo jRsJj^a 



Combining fl53|) . f lMjl . fl55|l imphes 

mF"(t)e(/i-5(t)) dpdxds < C [ [ [ [ F'{T)^{h-6{T)) dpdxdsdr 

and by Gronwall's lemma we deduce that 

Iff F'{t,s,x,p)^{h{s,x,p) -6{t)) dpdxds = 0, t eR+ 

and therefore 

snppF%t)c{{s,x,p) : h{s,x,p) < 6{t)} C {{s,x,p) : |p| < !^iNJ^ik!l|a;| < 

□ 

Once we have estimated the error between and /°, the asymptotic behaviour of 
f^{t,x,p) — f'^(t,t/e,x,p) as e \ follows by using the Sobolev inequality in Propo- 
sition 13.31 

Proof, (of Theorem 12.31) By Proposition 13.31 we have 

writ, ; ■) - f{t,t/e, OIIl^Cm^xms) = \\F%t, t/e, -, ■) - fit, t/e, -, OIU^CRBxRa) 

< C{m\F^(t)-f\t)\\ + \\TF^{t)\\) 



since Tf yt) = for any t G M_|_. Thanks to Proposition 16. II we know that 

\\F^it)-fm<C,{I)e, tG [0,/], e>0 
and we are done if we can find a constant C2{I) such that 

||rF^(t)|| <C2{I)e, te[0,I], e>Q 
since in that case we would obtain 

||r(t,-,-)-/°(t,t/£,-,-)llL2(M3xM3) <C(r)(Ci(/) + C2(/)K te[o,i], e>o. 
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Obviously, multiplying (|50|) by F*^ and integrating with respect to {s,x,p) we control 
the norm of F^{t, ■,-,■) uniformly in t G M+ and e > 

||F^(t)f = ||F^(0)|p = r||r||i.(^3,i,3). 

We intend to control the derivatives V(^t,x,p)F'^ as well, uniformly with respect to e > 0. 
The idea is to use the derivations commuting with T, introduced in Proposition 14.61 
Indeed, with the notation a(t,s,x,p) = {0, ^,eE(t,x) + Uc9{s) ^p) the equation (l50i) 
becomes 

dtF' + a ■ V(,,,,p)F^ + -TF' = 0. 

e 

Applying the operator d ■ V(s,x,p), i ^ {l;---;6} and taking into account that c* • 
^{s,x,p)TF^ = T{d ■ V(s,x,p)F^) one gets 

dtGt + a ■ V(^s,x,p)GI + [c\ a] ■ V^s,x,p)F' + -TGI = (56) 

where Gf = ■ V [s,x,p)F^ and [c*, a] are the Poisson brackets between the fields c', 
i G {!,..., 6} and a. Multiplying (15^ by and integrating with respect to {s^x^p) 
yield 

~fl I \Gl?<^pMs = - n [ GUc\a]-Vis,x,p)F' dpdxds 
^ Jo Jk^Jk'^^ Jo Jr^Jrs 

< \\Gm\\ ||[c\a]-V(,,.,p)F^(t)||, zG {!,... ,6} 

or equivalently 

l|G'^WII< 11^^(0)11+ ri|[c\a]-V(.,.,,)F^(r)l|dr, zg{1,...,6}. (57) 

Jo 

It is easily seen that for any i G {!,... ,6} the field [c*,a] has no component along s, 
since c* = = 0. Therefore [c\ a] G span{c^, c®} 

6 

[c\a\ = ^-fij{t,s,x,p) c?, % G {!,..., 6} 

i=i 

for some coefficients G Lii„^(]R+; L°°(R x x rs^), i,j g {!,..., 6}. Actually we 
have 

[c^ a] = ^^c^ + ec?,,E,c^ + (e^.^E^ - ^^^(.)) + e^.^^sc^ 
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[c\a] = J2ed.,E,c^^^ [c^ a] = --^^(.) c^ [c^ a] = -+^^(.) c^ [c^a] = -. 
^-^ ml ml m 



It follows that [c\a] ■ V [s,x,p)F^ = Yl^j=ilifi) ^'^'^ ( 157|) yields 



t 6 







We deduce that 



i=l j=l "^"^ i=l 

with 7(r) = maxjg|i_..._e} XlLi ll7ij('^)IU°° ^^'^ t)y Gronwall's lemma we obtain that 

{V(,,p)F^(t) : t G [0,/], £>0} 

remains bounded in L|^(Ms; L^(Ri^x]Rp)). By Proposition l6.2l we know that Ue>o,te[o,/]Supp F'^it) 
remains into a compact set of R^/TZ x R'^ x R'^ and clearly there is a constant C^{I) 
such that for any t G [0, /] 
P 



sup 

£>0 



V.F'{t) + {eE + uj,e{s)^p)-V:,F' < Cs{I){l + \\E{t)\\Lo.^^s)) (5^ 
m 



< C4(/)(l + 11^(0) |Uoo(K3) + \\dtE\\Lli[0,t];L^iR^))). 

It remains to estimate the time derivative dtF"^. As before we write 

dt{dtF') + a ■ Vis,x,p)idtF') + edtE ■ VpF' + ^r{dtF') = 

implying that 

liTMF'f = -e /7 / dtE{t)-VpF%t)dtF'{t)dpdxds 
1 Jo JrsJks 



< \\edtE{t)-VpF' 



We deduce that 



\\d,F%t)\\ < \\d,F%0)\\ + f \\ed,E{T)\\,^\\V,F%T)\\ dr. 

The family of time derivatives {dtF'^{t) : t G [0,/],£: > 0} remains bounded in 
L|(R,; L2(R3 X R3)) iff {\\dtF%0)\\ : e > 0} remains bounded in L|(R,; L\M.l x R^)). 
Notice that T-F^(O) = and therefore 

snp\\dtF%0)\\ = sup||-a-V(.,.,,)F^(0)|| (59) 

e>0 e>0 



sup 

e>0 



-a . V(,,.,,)r {x,p- ^eis) + ^^(0) ^x) 



< +00. 
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Combining fl58|) . f l59|l we deduce that there is a constant C2(/) such that 

(^t + £ ■ V,. + {eE{t) + u;,e{s) ^p) ■ V^) F%t)\\ < eC^il) 
saying that supg>o,te[o,/] < ^C'2(/)- □ 



sup 

£>o,te[o,/] 



7 Three dimensional setting 

In this section we study the particle dynamics under fast oscillating three dimensional 
magnetic fields 

B'{t,x) = e{t/e)B{x)b{x), div^.(56) = 

for some scalar positive function B{x) and some field of unitary vectors b{x) G M^. 
The analysis is completely analogous to that for fast oscillating homogeneous magnetic 
fields previously discussed. Therefore we only focus on the formal derivation of the limit 
model. By Gauss's magnetic law div^.i?^ = we can write Bb = curl^^A, div^.^ = 
and by Faraday's law dtB'^ + curl^E'^ = we deduce that the rotational part, curia-?/', 
of the electric field = — Vx0 + curia;?/' is given by 

curla.V' = -^9'{t/e)A{x). 



The Vlasov equation becomes, with the notations E = —'Vx4>, 0Jc{x) 



eB{x) 



m 



dtf + - ■ V.r + {eE{t, x) - ^l^^^A{x) + u,{x)e{t/e)p A b{x)\ ■ V.f = 0. (60) 
We prescribe the initial distribution 

r(0,a;,p) = r(a;,p), (a;,p) G x Ml (61) 
Plugging the Hilbert expansion ([TT]) into (IBU|1 yields 

dsf-ee\s)A{x)-Vpf = Q (62) 

at the lowest order and 

dtf + - ■ Va./° + {eE{t, x) + u,{x)e{s)p A b{x)) ■ + dj^ - e 9 '(s) A(x) ■ V^/^ = 
m 

(63) 
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at the next order e^. As before, the point is how to ehminate the first order correction 
appearing in (jHS]), based on the constraint f p^ . We introduce the operator 

Tiu = div(s,p){M(l, -e 9 '{s)A{x))} 

with domain 

DiTi) = {u e LliRs-, L^RlxRl)) : div(,,p) {u (1, -e ^ '(s)A(x))} G L|(M,; ^^(M^ xMJ))}. 

The characteristics {S,X,P){t; s,x,p) of the first order differential operator 7i are 
given by 

S{t] s,x,p) = s + T, X{t] s,x,p) = X, P{t] s,x,p) = p + e{9{s) — 9{s + t))A{x). (64) 

Notice that a complete family of functional independent invariants is given by {x,p + 
e9{s)A{x)} and therefore the constraint (162 p becomes 

3g' = g%t,x,q) : f{t, s,x,p) = g%t,x,q = p + e9{s)A{x)). 

In particular f^{t,t/e, x,p) is fast oscillating through the periodic profile 9{s = t/e) and 
therefore we expect that f^{t, s,x,p) is the two-scale limit of f^{t,x,p) when e \ 0. 
The average operator (■)^ along the characteristic flow f l64|) is given by 

1 /-^ 

{u)i{s,x,p) = f u{S{t;s,x,p),X{t;s,x,p),P{t;s,x,p)) dr (65) 

u{s + T,x,p + e{9{s) - 9{s + t))A{x)) dr 

M (r, X, p + e 6'(s)yl(x) - e 6'(r) A(x)) dr 

for any function u G L^(R<j; -^^^(M^ x M^)). The dynamics for /° is obtained by elim- 
inating in f l63|) taking into account that the functions in the range of 7i are zero 
average. We have 

dtf + - ■ V,f + {eE{t, x) + uj,{x)9{s)p A h{x)) ■ V^f = -TJ' G Range Ti = ker (■), 
m 

and therefore (l63l) is equivalent to 

(dtf + - ■ V./° + {eE{t, x) + u,{x)9{s)p A b{x)) ■ V^f) = 0. 
\ m / 1 

We need to average the derivatives with respect to (t, x,p) of the density /°, under the 

constraint (162|) . Clearly we have {dtf)i = dt{f)i = dtp. 
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Lemma 7.1 Assume that f{s,x,p) = g {x,q = p + e 9{s)A{x)) is smooth. Then we 
have 

■ Vj) = ^ ~ ' "^^^^ ■ V.g + -d^Am q - e {6') A) • V,g 
\m / 1 m m 

= p + . V./ + ^d^A {m - e)p + e{2e (0) -e'- {e'))A{x)) ■ v,/ 

and 

{{eE{x) +uj,e pAb) -Vpf)^ = [eE + ujc{{e)qAb-e{e^)A{x) Ab)] -Vqg 

= [eE + ujc{e)pAb + uj^e{e {6} - {e^))A{x) A b] ■ VJ. 

Proof. We have 

V,/ = V,g + e 9{s) 'd.AVgg, Vpf = V^g. 

Since V {x,q)g are constant along the flow flM|l we can write 

(- ■ V./\ = Mi . v,(? + - (p . 'd:,AV,g. 
\m / 1 m m 

It is easily seen by the definition of the average operator (■)^ that 

(p)^ =p + e e{s)A{x) - e (9) A{x) 

and 

(p 6)^ = {p + e e{s)A{x)) (6) - e {6^) A{x) 

implying that 

■ V./) = ^ ~ ' ^^""^ ■ V.g + -d.Am q - e {9^) A{x)) ■ V,g. 
\m / 1 m m 

Similarly one gets 

((eE(x) +UcepAb)- Vp/)i = [eE + qAb-e (6^) A{x) A b)] ■ V^g. 

□ 

Combining the previous computations and using the identities 

dr^A A + AA cnil^A = ^d^A A, d^A q + q A curl^A = ^d^A q 
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yield the transport equation in the space phase (x, q) 

^^^0^ g-^j^MM. v,/+ (eE + - (9) %A q - - {9') %A a) -V.g' = (66) 

where s, = g^it, x,q = p + e9{s)A{x))^ since 7i/°(t) = for any t G M_|_. The 
transport equation in the phase space becomes 



p + e{9{s) - {9))A{x 



m 



+ L^^l^^^lll^^^lll^.^jo^ eE + ^m'd.A-9d.A)p 



m 



VX = 0. (67) 



+ -{9 (9) - {9')) 'd,AA + -{{9) - 9) 9 d.,AA 
m m 

We supplement these transport equations by the initial conditions 

g\0,x,q) = rix,q-e9iO)Aix)) (68) 

f (0,s,x,p) = r(x,p + e (0(s) -0(O))A(x)). (69) 

Following the lines in Sections El IH] we can prove weak and strong convergence results, 
which justify the Hilbert expansion in f[T7|) . In the weak framework we obtain 

Theorem 7.1 Assume that E e Ll^(R+; L°° (R^)) , A e Ll^{R+]W^'°°{R^)f , G 
L2(M3 X R^). For anye>0 let f G L°^(M+; L'^{R^ x R^)) be a weak solution of §0!), 
(l6Ti) . Then there is a sequence En \ such that {f'^")n two-scale converges towards a 
weak solution of (jEZD; dSHD- 
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